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Saccadic  eye  movement*  are  among  the  fastest  voluntary  muscle  movements 
the  human  body  is  capable  of  producing  and  are  characterized  by  a  rapid  shift 
of  gaze  from  one  point  of  fixation  to  another.  Although  the  purpose  for  such 
an  eye  movement  is  obvious,  that  is,  to  quickly  redirect  the  eyeball  to  the 
target,  the  neuronal  control  strategy  is  not.  For  instance,  does  the  word 
"quickly"  in  the  previous  sentence  imply  the  most  rapid  movement  possible,  or 
siaqily  a  fast  movement  as  opposed  to  a  slow  movement?  This  investigation  uses 
system  identification  techniques  to  estimate  muscle  forces  during  horizontal 
saccadic  eye  movements  in  order  to  better  understand  the  neuronal  control 
strategy. 

Until  recently,  sradels  of  the  saccadic  eye  movement  system  involved  a 
ballistic  or  preprogrammed  control  to  desired  eye  position  based  on  retinal 
error  alone  (1-4).  Today,  an  increasing  number  of  authors  are  putting  forth 
the  idea  that  goal-directed  saccades  are  controlled  by  a  local  feedbaok  loop 
that  continuously  drives  the  eye  to  the  desired  position.  This  hypothesis, 
first  presented  by  Vossius  (5)  in  1960,  did  not  start  to  gain  acceptance  until 
1975  when  Robinson  (6)  reexamined  it.  Robinson  suggested  that  sacoades 
originate  from  neural  commands,  which  encode  saocade  velocity  and  duration,  to 
the  pulse  generator  that  specifies  the  desired  position  of  the  eye  rather  than 
the  preprogrammed  distance  the  eye  must  be  moved.  The  image  of  the  internal 
representation  of  the  current  eye  position  is  subtracted  from  the  internal 
representation  of  the  desired  position,  creating  an  error  signal  that 
completes  the  local  feedback  loop  to  generate  the  neural  control  signal.  This 
neural  pulse  continuously  drives  the  eye  until  the  internal  representation  of 
the  error  signal  is  zero. 

One  of  the  basic  tasks  of  this  project  is  to  investigate  whether  an 
external,  continuous-time  sensor  feedback  control  mechanism  operates  during  a 
saccadic  eye  movement.  Studies  by  other  investigators  indicate  that  the 
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profile  of  i  seceede  can  be  altered  up  to  50  ma  after  target  aoveaent,  and 
that  visual  information  is  processed  continuously  during  this  interval  (7,8). 
Their  results  indicate  that  50  as  after  target  aoveaent,  the  saccade  profile 
cannot  be  altered  by  altering  the  target  location.  This  ezperiaent  atteapts 
to  verify  their  results  by  artificially  aaintaining  a  constant  error 
throughout  the  saccade  by  moving  the  target  the  saae  distance  the  eyeball 
aoves.  Collevijn  and  Van  der  Mark  (9),  in  their  study  of  the  slow  phase  of 
optokinetic  nystagaus,  used  this  aethod  of  artificially  opening  the  feedback 
loop. 

A  second  task  in  this  project  involves  the  investigation  of  a  horizontal 
saccadic  eye  aodel  presented  by  Enderle  (10).  Estiaated  values  of  the  poles 
of  the  transfer  function  shoved  an  increasing  trend  as  target  displaceaent 
increased.  A  nuaber  of  reasons  aay  account  for  this  increase,  a  significantly 
nonlinear  aodel  and/or  an  inaccurate  awdel,  or  poor  initial  estimates  for  the 
paraaeter  values.  The  aodel  for  horizontal  saccadic  eye  movements  must  be 
accurate  since  it  will  be  used  in  sensitively  assessing  whether  an  external 
feedback  control  aechanisa  operates  during  saccade s. 


MHO.  OF  OCDLOHOTOK  PLANT 

In  developing  Enderle* s  aodel  for  horizontal  saccadic  eye  aoveawnts, 
equation  2  assumed  that  the  tine  derivative  of  the  antagonist  tension  was 
negligibly  saall  (10).  This  assumption  was  a  aajor  factor  in  the  resulting 
functional  relationship  between  target  displaceaent  and  poles  of  the  transfer 
function.  The  initial  paraaeter  estimates,  which  will  be  discussed  in  the 
next  section,  also  affect  this  relationship.  Using  Laplace  variable  analysis 
about  the  operating  point  or  initial  eye  position,  it  is  possible  to  derive 
the  following  new  and  aore  accurate  aodel  of  horizontal  saccadic  eye  aoveaents 
appropriate  for  the  systea  identification  technique  (11). 
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This  aodel  aodifies  and  corrects  the  linear  hoaeoaorphic  aodel  by  Bahill  et 
al.  (4).  Iaplementing  this  aodel  in  the  paraaeter  estimation  routine 
resulted  in  an  excellent  aateh  between  aodel  prediction  and  the  data. 


Figure  1  presents  s  block  diagram  of  the  saccadic  eye  movement  system 
which  functions  as  one  element  within  the  overall  system  controlling  all 
oculomotor  movements.  The  saccadic  eye  movement  system  consists  of  three 
components:  (1)  the  controller  (i.e.,  the  agonist  and  antagonist  active  state 
tension  generators),  (2)  the  oculomotor  plant,  and  (3)  the  feedback  element 
H.  The  oculomotor  plant,  indicated  by  the  section  within  the  dashed  lines  in 
Figure  1,  is  based  on  the  Laplace  transform  of  equation  1.  Note  that  under  an 
assumed  ballistic  control  of  saccadic  eye  movements,  time  delays  in  the  feed¬ 
back  element  cause  this  system  to  operate  in  an  open-loop  mode.  If,  in  fact, 
an  external  continuous-time  feedback  control  mechanism  exists,  then  opening 
the  feedback  loop  causes  the  system  to  overreact  because  of  the  constant  error 
present  at  the  input. 


Figure  1.  A  block  diagram  of  the  oculomotor  plant.  The  section  of  the 
disgram  within  the  dashed  lines  is  the  modified  linear 
homeomorphic  model.  The  feedback  element  H  is  unity. 

Since  the  derivation  of  the  modified  linear  homeomorphic  model  requires  no 
assumptions  concerning  the  saccadic  active  state  tensions,  any  type  of  control 
signal  may  be  used.  For  example,  investigators  may  use  this  model  when 
examining  alternatives  to  the  pulse-step  muscle  force.  If  a  continuous-time 
external  feedback  control  mechanism  exists,  then  a  new  waveform  could  be 


postulated  for  the  active  state  tension.  This  waveform  coaid  then  be  revised 
based  on  the  predicted  response  of  the  modified  linear  homeomorphic  model. 

The  mechanical  components  of  the  model  are  connected  so  that  the  active 
state  tension  rate  of  change,  together  with  the  actual  active  state  tensions, 
drives  the  eyeball  to  its  final  position.  Consider  the  low-pass  filtered 
pulse-step  active  state  tension  signals  presented  in  Figure  2.  During  the 
initial  pulse  phase  of  the  trajectory,  the  agonist  active  state  tension  drives 
the  eyeball  and  the  antagonist  active  state  tension  restrains  the  eyeball. 
Since  the  antagonist  active  state  tension  falls  to  zero  during  this 
interval,  then  the  antagonist  active  state  tension  rate  of  change,  F^^,,  is 
negative  and,  as  indicated  by  equation  1,  acts  to  drive  the  eyeball  along 
with  the  agonist  force,  F ^  ,  and  the  agonist  rate  of  change  force,  F^,.  As  a 
result,  a  greater  force  propels  the  eyeball  than  if  only  the  agonist  and 
antagonist  active  state  tensions  are  nsed,  as  shown  in  the  plant  input  curves 
of  Figure  2.  This  result  agrees  with  the  tension  data  recorded  for  horizontal 
saccadic  eye  movements  (12,13).  After  the  pulse  phase  of  the  trajectory,  the 
agonist  active  state  tension  falls  to  a  steady  state  value,  resulting  in  a 
negative  agonist  active  state  tension  rate  of  change,  fag-  Thus,  according  to 
equation  1,  the  agonist  rate  of  change  force,  FAG'  acts  to  restrain  the 
eyeball  with  the  antagonist  force,  F^^.,  and  the  antagonist  rate  of  change 
force,  F.„_.  As  shown  for  the  plant  input  curves  in  Figure  2,  the  resultant 
forces  act  like  a  dynamic  brake  which  is  reduced  to  zero  as  the  eyeball 
reaches  its  destination. 


SYSTEM  IDENTIFICATION  TECHNIQUE 

Parameter  estimates  for  the  oculomotor  system  model  presented  in  the 
previous  section  are  found  using  the  system  identification  technique  as 
described  by  Enderle  (10,14).  The  system  identification  technique  is  a 
frequency  response  method.  Ideally,  the  transfer  function  is  equal  to  the 
ratio  of  the  Fourier  transform  of  the  output  to  the  system  input.  For  the 
oeulomotor  system,  the  transfer  function  is  calculated  from  the  fast  eye 
response  to  a  step  target  displacement.  Unequal  time  delays  and  eye 
displacements  from  saccade-to-saccade  for  the  same  target  displacement,  make 
it  impossible  to  use  averaging  techniques  to  reduce  the  effects  of  measurement 
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Figure  2.  Transformation  of  the  active  state  tensions  by  the  numerator 
terms  of  the  oculomotor  plant.  One  diagram  shows  the  assumed 
low-pass  filtered  pulse-step  waveform  of  the  agonist  and 
antagonist  active  state  tensions  initiated  at  100  ms.  The 
antagonist  pulse  circumscribes  the  agonist  pulse  by  3  ms  on  each 
end.  Parameter  values  are  from  Bahill  et  al.  (4)  for  simulating 
a  10  horixontal  saccadic  eye  movement.  The  other  diagram 
ahows  a  partial  and  full  transformation  of  the  active  state 
tensions  by  the  numerator  terms  of  the  oculomotor  plant. 


noise.  Fortunately,  the  measurement  noise  is  small  relative  to  the  input  and 
output  signals,  so  the  second  method  described  by  Otnes  and  Enochson  is 
used  (15).  The  fast  eye  response  measurements  are  first  filtered  using  a 
Butterworth  low— pass  filter  with  a  half-power  point  at  125  Hz.  Transforming 
the  filtered  measurements  directly  by  the  Fast  Fourier  algorithm  resulted  in 
distortion  due  to  truncation,  since  the  signal  did  not  go  to  zero  at  steady 
state.  This  problem  is  circumvented  by  subtracting  the  steady  state  value 
from  each  sasiple,  passing  this  signal  through  a  Kaiser  window,  packing  with 
zeros,  and  Fast  Fourier  transforming  the  modified  signal.  The  Fourier 
transform  of  the  fast  eye  response  is  equal  to  the  Fast  Fourier  transform  of 
the  modified  signal  plus  the  Fourier  transform  of  the  unit  step  with  amplitude 
equal  to  the  steady  state  value.  The  input  signal  is  the  Fourier  transform  of 
the  unit  step  function  with  amplitude  equal  to  the  steady  state  value. 
Parameter  estimates  for  the  oculomotor  model  are  calculated  using  the 
conjugate  gradient  search  program  similar  to  Seidel's,  which  minimizes  the 
integral  of  the  absolute  value  of  the  error  squared  between  the  model  and  the 
data  (16). 

Data  were  collected  from  subjects  seated  before  a  target  display  of  seven 
small  red  light-emitting  diodes  (LED),  each  separated  by  5°.  The  subject's 
head  was  restrained  by  a  bite  bar.  The  subject  was  instructed  to  follow  the 
"jumping"  target  which  moved  from  the  center  position  to  any  of  the  other  LEDs 
and  returned  to  the  center  position.  Target  position  order  and  time  interval 
between  displacements  were  randomized.  Data  were  only  recorded  for  the  initial 
displacement  from  the  center  position. 

Horizontal  eye  movements  were  recorded  from  each  eye  using  an  infrared 
signal  reflected  from  the  anterior  surface  of  the  cornea-scleral  interface 
with  instrumentation  described  by  Engelken  et  al.  (17).  Signals  for  both  eyes 
tracking  were  digitized  using  the  analog/digital  converter  of  the  DECLAB  PDP 
11/34  computer  and  stored  in  disk  memory.  These  signals  were  sampled  at  a 
rate  of  1000  samples  per  second  for  one-half  second. 

Figure  1  presents  a  block  diagram  of  the  oculomotor  system,  which  is 
suitable  for  the  case  of  a  ballistic  controller,  or  as  indicated,  with  an 
external  continuous-time  feedback  control  configuration.  Under  the  assumption 
of  ballistic  control,  the  system  operates  in  an  open-loop  mode  due  to  time 
delays  in  the  feedback  element.  Under  the  assumption  of  an  external 
continuous-time  controller  during  experiments  in  which  the  feedback  loop  is 


artificially  opened,  the  system  also  operates  in  an  open-loop  mode.  Thus, 
the  analysis  of  the  data  is  performed  with  an  open-loop  system  regardless  of 
the  type  of  control  assumed  operating.  However,  the  signals  describing  the 
agonist  and  antagonist  force  generators  are  quite  different,  dependent  on  the 
type  of  controller. 

The  following  low-pass  filtered  pulse-step  waveforms  describe  the  agonist 
and  antagonist  active  state  tensions  under  an  assumed  ballistic  control 
mechanism: 
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initial  magnitude  of  the  agonist  active  state  tension 

pulse  magnitude  of  the  agonist  active  state  tension 

step  magnitude  of  the  agonist  active  state  tension 

initial  magnitude  of  the  antagonist  active  state  tension 

step  magnitude  of  the  antagonist  active  state  tension 

latent  period,  the  time  interval  between  the  target  movement  and 

the  start  of  the  eye  movement 

duration  of  the  agonist  pulse  active  state  tension 
activation  time  constant 
deactivation  time  constant. 


The  steady  state  antagonist  active  state  tension  terms,  and  are 

removed  from  the  analysis  since 


and 
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according  to  equation  1.  If  either  active  state  tension  is  less  than  zero,  it 
is  set  equal  to  zero. 

Physiological  data  was  used  to  estimate  all  of  the  initial  parameter 
estimates  (4,13,18).  Specifically 


V 


I  -  66.4  N/a 
I^T  *  32.  N/m 

KSE  =  125  N/" 

=  3.5  N  s/n 

BANT  =  1’66  N  s/“ 

B  =  3.1  N  s/m 

J  -  2.2xl0"3  N  s2/i 


(Robinson) 
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The  value  of  K  is  determined  from  the  steady  state  agonist  and  antagonist 


muscle  tensions 


where 


T  -  T  =  tr 
AG  ANT 


T^  **  17  +  ©  g  tension 


T^j,  =  17  -  0.3  ©  g  tension 

from  Figure  1  of  Robinson  et  al.  (18).  The  value  of  the  muscle  viscosity 
terms  reported  by  Bahill  et  al.  are  modified  according  to  their  incorrectly 
derived  result  (4). 

bag  *  kst  bag  /  kse 

and 

BANT  =  KST  BANT  /rSE 

Substituting  these  values  into  the  expressions  for  6,  Pq,  P^,  Pj  and  P^  found 
in  Enderle  et  al.  (11),  yields 


=  5.0938109  x  10' 

=  1.548  x  103 

=  3.4455  x  105 

=  1.9724  x  107 

=  2.2631  x  108 


The  initial  estimates  for  the  four  real  eigenvalues  as  determined  from  the 
characteristic  equations  for  this  system  are  -15.27,  -66.1,  -173.36,  and 
-1293.28.  Now,  the  steady  state  active  state  tensions,  F^q,  ^AGS'  are 
computed  from 
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mm 
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(3) 


PAG  "  Tag  +  ^(x  -  xpi) 

SE 


derived  froa  Figure  1  in  Enderle  et  al.  (11).  Substituting  equation  2 
equation  3  with  xpj  equal  to  3.14  aa  (assumed  from  Robinson  (18))  yields 


AG 


.1089  +  .01850  N 


into 


The  initial  estiaate  of  the  pulse  magnitude  of  the  active  state  tension,  F^op, 

is  0.9806  (a  maximum  innervation  signal  of  100  g)  (13)  . 

The  initial  estiaate  of  the  duration  of  the  agonist  pulse  active  state 

tension  is  estimated  directly  from  the  data.  The  maximum  velocity  is  computed 

using  a  two-point  central  difference  derivate  method  (19).  The  initial 

estimate  of  t  is  then  the  time  interval  between  the  start  of  the  saccade  and 
P 

the  time  at  maximum  velocity. 

The  initial  estimate  of  the  deactivation  time  constant  is  also 
estimated  directly  from  the  data  as  the  time  interval  between  maximum  velocity 
and  the  end  of  the  saccade  divided  by  4.  The  activation  time  constant  equals 
2  ms  (20). 


Great  care  was  exercised  in  evaluating  initial  parameter  estimates, 
since  large  differences  from  the  true  values  conld  cause  the  estimation 
routine  to  converge  to  suboptimal  and  nonphysiologically  consistent  results. 
The  precision  of  the  parameter  estimation  routine  results  presented  in  Figures 
3  and  4  for  a  10°  target  movement  are  typical  for  the  three  subjects  tested. 
There  was  no  significant  trend  in  the  values  of  the  eigenvalues  with  target 
displacement  and  the  estimates  correlated  well  with  their  physiological  data 
derived  values.  Table  1  lists  the  final  estimated  values  for  the  eigenvalues 
for  target  displacements  of  5,  10,  and  15°  for  one  of  the  subjects. 


EXTERNAL  FEEDBACK  CONTROL  MECHANISM 
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Experiments  were  conducted  in  which  the  hypothesized  continuous-time 
external  feedback  control  mechanism  was  artificially  opened.  Data  were 
collected  on  three  subjects  in  two  stages.  First  stage  data  consisted  of 
natural  saccades,  elicited  from  target  displacements  of  -15  to  15°.  The 
second  stage  consisted  of  data  collected  by  artificially  maintaining  a 
constant  error  of  5°  throughout  the  saccade  by  moving  the  target  (under 
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TABLE  1.  FINAL  PARAMETER  ESTIMATES  FOR  THE  FOUR  REAL  EIGENVALUES  FOR  TARGET 
DISPLACEMENTS  OF  5,  10,  AND  15°  LEFT  OF  CENTER  POSITION 


TARGET  DISPLACEMENTS 


5° 

10° 

15° 

LEFT  BIB 

Eigenvalue  1 

-14.4 

-13.9 

-13.7 

Eigenvalue  2 

-83.8 

-84.4 

-93.6 

-190.9 

-174.6 

-157.4 

Eigenvalue  4 

-1298.7 

-1283.2 

-1259.4 

RIGHT  BXB 

Eigenvalue  1 

-14.8 

-14.4 

-12.35 

Eigenvalue  2 

-80.8 

-84.9 

-88.2 

Eigenvalue  3 

-183.8 

-165.4 

-160.1 

Eigenvalue  4 

-1302.2 

-1261.6 

-1223.7 

computer  control)  the  sane  distance  the  eyeball  moved.  This  type  of  movement 
will  be  referred  to  as  the  modified  target  trajectory.  The  experiment  started 
by  the  target  mowing  +5°.  The  subject  responded  by  following  the  jumping 
target.  Instantaneously,  as  the  eyeball  moved,  the  target  moved  the  same 
distance.  Significant  differences  in  saccade  profile  between  the  two  stages 
would  indicate  the  possible  existence  of  an  external  feedback  mechanism.  Data 
was  first  analyzed  using  the  two-point  central  difference  method,  and  then 
using  the  system  identification  technique. 

The  results  of  the  analysis  of  the  data  on  the  three  subjects  using  the 
two-point  central  difference  method  are  contained  in  the  Appendix.  Analysis 
results  did  not  support  the  postulated  continuous-time  external  feedback 
control  mechanism.  However,  the  time  interval  from  the  start  of  the  saccade 
to  the  time  at  peak  velocity  for  subjects  2  and  3  did  show  marked  variation 
ranging  from  21  ms  to  8  ms.  Most  of  the  values  were  clustered  about  evenly 
into  two  groupings  of  16  ms  and  10  ms.  These  time  interval  results  are 
plotted  in  Figure  5.  It  is  interesting  to  note  that  the  5°  movements  elicited 
under  natural  conditions  were  always  in  the  larger  time  interval  grouping  of 
16  ms  for  each  of  the  three  subjects.  It  is  speculated  that  if  there  were  a 
continuous-time  external  feedback  control  mechanism,  then  both  the  time 


interval  and  the  final  eye  displacement  would  have  been  noch  larger  than 
under  normal  conditions. 
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Figure  5.  The  time  interval  from  the  start  of  the  saecade  to  the  time 
at  peak  velocity  vs  absolute  value  of  displacement  for  sub¬ 
ject  2. 


The  results  of  the  analysis  of  the  data  on  the  three  subjects  using  the 
system  identification  technique  confirmed  the  previous  findings. 
Additionally,  estimates  of  the  pulse  magnitude  of  the  agonist  active  state 
tension  F^p  were  clustered  about  evenly  into  two  groupings,  corresponding  to 
the  previous  groupings,  0.95  N  and  0.85  N.  These  results  are  plotted  in 
Figure  6.  The  estimated  pulse  magnitude  for  the  agonist  active  state  tension 
as  a  function  of  displacement  for  the  natural  saccades  is  illustrated  in 
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Figure  6.  The  estimated  values  of  the  pulse  magnitude  of  the  agonist  active 
state  tension  vs  absolute  value  of  displacement  for  subject  3 
during  the  stage  2  experiment. 
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Figure  7.  The  estimated  values  of  the  pulse  magnitude  of  the  agonist  active 
state  tension  vs  absolute  value  of  displacement  for  subject  3 
during  the  stage  1  experiment. 


Figure  7.  One  prononnced  feature  evident  from  these  two  graphs  is  that  there 
does  not  appear  to  be  a  strong  relationship  between  pnlse  magnitude  and 
displacement.  In  faot  the  variation  in  F^,p  seen  in  Figure  6  is  approximately 
the  same  size  variation  seen  in  Figure  7.  Apparently  the  magnitude  of  the 
agonist  pulse  is  maximum,  regardless  of  the  amplitude  of  the  saccade,  and  only 
the  duration  of  the  agonist  pnlse  affects  the  size  of  the  saccade.  Under 
these  conditions,  the  eyeball  is  driven  to  its  destination  in  minimum  time  for 
saccades  of  all  sizes  (21). 


CONCLUSION 

The  conclusions  obtained  during  the  performance  of  this  project 
indicate  that  the  saccadic  eye  movement  system  operates  under  a  time  optimal 
control  strategy  without  continuous- time  external  feedback  sensory 
information.  The  implication  that  the  eyeball  is  directed  to  its  destination 
in  minimum  time  for  saccades  of  all  sizes  certainly  deserves  additional 
attention,  considering  the  small  number  of  subjects  tested.  This  new 
hypothesis  certainly  should  be  investigated  further  in  light  of  the  current 
theory  on  the  neuronal  control  of  saccadic  eye  movements.  Particular 
attention  should  be  paid  to  the  large  variation  seen  with  saccades  of  the  same 
size.  Furthermore,  the  importance  of  the  second  phase  of  the  saccade,  that 
is,  when  the  agonist  pulse  is  falling  to  its  steady  state  value,  must  be 
understood  in  relation  to  the  time  optimal  approach. 
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ANALYSIS  RESULTS 


The  two-point  central  difference  analysis  resalts  on  three  subjects  for 
natural  and  aodified  target  trajectory  saocade  data  are  shown  in  Tables  A-l 
through  A-6. 

List  of  Symbols 

Saeeade  -  distance  the  target  moved 

Left  and  Right  -  the  distance  the  left  and  right  eye  moved  in 
Displaeeaemt  response  to  the  target  displacement  (degrees) 

IC  -  initial  position  of  the  eyeball  (degrees) 

TO  -  time  from  the  start  of  the  saccade  to  the  time  at  peak 
velocity  (ns) 

DOR  -  duration  of  the  saccade  (ms) 

TP  -  time  at  which  the  saccade  started  (ms) 

PV  -  peak  velocity  of  the  saccadic  eye  movement  (degrees/second) 

DI8T  -  eye  position  at  the  time  of  peak  velocity  (degrees) 


TABLE  A-l .  MODIFIED  TARGET  TRAJECTORY  DATA 
FOR  THE  LEFT  EYE  OF  SUBJECT  1 


